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1 .  Introduction 

We  consider  a  2m  factorial  experiment  under  a  completely  randomized 

design.  The  runs  are  denoted  by  t^  ■  (t^ . tiu,***,tinJ  w^ere  c^u  *-s 

the  level  of  the  factor  u  in  the  run  i,  t^u  •  0,1;  u  -  1 . . 

i  *  Let  T  (nxm)  be  a  matrix  whose  ith  row  is  t^,  i  *  l,...,n. 

The  matrix  T  is  called  the  array  (or  the  design)  for  the  experiment.  We 
take  r^(>  1)  observations  for  the  ith  (i  *  l,...,n)  run  of  the  experi¬ 
ment.  The  experiment  just  described  is  very  common  in  quality  control 
studies  (see  Taguchi  and  Wu  1985)  where  a  major  goal  is  to  evaluate  the 
sensitivity  of  the  manufacturing  process  to  noise.  The  runs  t] , 
i  *  l,...,n,  in  T  are  called  the  level  combinations  of  m  control 
factors.  An  important  problem  is  to  find  an  optimum  combination  of 
levels  of  control  factors  in  view  of  reducing  the  process  variability  due 
to  noise  factors.  A  list  of  noise  factors  likely  to  affect  the  process 
variability  is  first  made.  Various  level  combinations  of  noise  factors 
that  provide  a  good  representation  of  noise  are  then  selected.  For  the 
ith  run  _t^,  we  collect  r^  observations  corresponding  to  r^  level  combina¬ 
tions  of  noise  factors.  The  r^  observations  for  a  run  in  T  are  called  r^ 
replicated  observations.  The  variability  in  the  r^  replicated  observa¬ 
tions  for  the  ith  run  _tj  is  attributed  to  the  process  variability  due  to 
noise  in  the  experiment. 

Let  y^j  be  the  jth  observation  on  the  ith  run  t_|  in  T, 

j  -  l,...,rif  i  -  1 , . . .  ,n.  We  denote  jr  -  (yn,"',ylr  ’ '  *  * ;ynl ’ *  * *ynr  * 

1  n 

There  are  N  *  ri+*,,+rn  observations  in  The  model  for  the  experiment 


is 
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E(^)  -  XB,  (1) 

V<£>  “  S,  (2) 

where  X(Nxp)  is  a  known  matrix  that  depends  on  the  array  T,  J3(pxl)  is  a 
vector  of  factorial  effects  considered  in  the  experiment.  The  matrix  T. 
is  an  unknown  diagonal  matrix  where  the  diagonal  elements  for  all  obser¬ 
vations  on  *  (tu,...,tin)  are  equal  to  o^t^,...,^  )  «  o2(j^). 

The  parameters  o  (t^,...,t^u  *  )  are  called  the  dispersions  of 

the  factor  u(u  -  l,...,m)  at  level  1.  The  parameters  o2(til , . . . ,tiu  - 
are  called  the  dispersions  of  the  factor  u(u  ■  l,...,m)  at 
level  0.  We  are  interested  in  comparing  the  dispersions  at  level  1  with 
the  dispersions  at  level  0  for  every  factor.  For  this  purpose,  we  define 
the  dispersion  effect  of  the  factor  u(u  =  l,...,m)  as 
rc(t11,...,tiu  -  . tiu  - 

-£d(ttl,...,tlu  -  0»"*»tim^°2(ttl'*,,’tiu  “  0 . ciiJ»  (3) 

where  c( t^ ^ , . . . ,t^u  -  1 , . • . , t )  and  d(t^ ^ , . . . ,t^u  *  0,...,tlffl)  are  known 

constants  with  Ecf t ,t,  -  l,...,t,  )  -  £d(t.. . .  *  0,...,t,  ) 

il  iu  im  v  ll  iu  im' 

*  1.  The  dispersion  effect  is  a  contrast  of  the  dispersions  at  level  1 
and  level  0.  If  the  dispersion  effect  is  greater  than  0,  we  then  prefer 
level  0  over  level  1  of  the  factor.  If  the  dispersion  effect  is  less 
than  0,  we  then  prefer  level  1  over  level  0.  If  the  dispersion  effect 
is  zero  then  levels  1  and  0  are  equally  preferable.  Notice  that  the 
preference  of  one  level  over  the  other  level  depends  on  the  dispersion 
effect  considered  or  in  other  words,  on  the  values  of  c's  and  d’s.  The 
linear  functions  of  dispersion  effects  of  the  form  (3)  are  used  in 
measuring  how  the  dispersion  effect  of  one  factor  varies  at  different 
levels  of  the  other  factors  [see  Box  and  Meyer  (1986),  Srivastava 


(1987)].  In  this  paper  we  consider  the  problem  of  estimating  meaningful 
dispersion  effects  of  the  factors.  We  also  characterize  arrays  so  that 
the  estimation  of  meaningful  dispersion  effects  is  possible.  Although  in 
(3)  we  use  the  contrast  of  level  1  and  level  0  dispersions  for  the 
dispersion  effect,  one  may  also  take  the  ratio  of  two  dispersions 
considered  (see  Box  and  Meyer  (1986)).  The  use  of  any  of  these 
definitions  will  however  give  the  same  conclusion  since  each  of  them 
compares  the  dispersions  at  level  l  and  level  0.  We  assume  throughout 
the  paper  that  for  the  fitted  model  to  the  data  there  is  no  significant 
lack  of  fit. 

The  dispersion  model  and  dispersion  effects  are  considered  in  the 
work  of  Taguchi  (see  Taguchl  and  Wu  1985),  Box  (see  Box  and  Meyer  1986), 
Kacker  and  Shoemaker  (1986),  Phadke  (1986),  Ghosh  (1987),  Srivastava 
(1987),  Nair  and  Pregibon  (1988).  This  paper  makes  a  further 
contribution  to  this  area  of  research. 

2.  Estimation  of  Dispersion  Effects. 

We  choose  a  factor  out  of  m  factors  and  present  estimators  of  some 
meaningful  dispersion  effects  of  the  chosen  factor.  For  clarity,  we  do 
not  introduce  any  notation  for  the  chosen  factor.  We  define  the 
following  indicator  variable  for  i  *  l,...,n,  so  that  it  is  possible  to 
identify  the  runs  at  which  the  factor  appears  at  level  1  and  level  0, 
respectively: 


6 


i 


1  if  the  level  of  the  factor  in  the  ith  run  is  1, 
, 0  if  the  level  of  the  factor  in  the  ith  run  is  0. 


A  meaningful  dispersion  effect  of  the  form  (3)  is 
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Z  6  a2(t’)  Z  (l-Sja2^’) 

a2(l)  -a2(0)  -ill- _  _ (4) 

n  n 

Z  5.  Z  (1-6  ) 

1-1  1  1-1 

The  dispersion  effect  (4)  is  the  difference  between  the  average  dis¬ 

persions  at  level  1  and  level  0  of  the  chosen  factor.  Let  y^  be  the 
mean  of  all  observations  on  the  ith  run,  i  -  l,...,n.  The  sum  of  squares 


of  pure  error  is 


n  ri 


n  ri 


n  ri 


SSPE  »  £  Z  (y  -  y. )2  =  Z  £6  (y  -;.)2+  Z  Z  ( 1  -  6  ) (y  -  y  )' ‘ 
i-11-1  1J  1  i-11  =1  1  1  i  —  1 1  —  1  1  J 


We  now  write 


n  ri 

r  z  6  (y  -  y  )2 

*}«>  — - - 

L  n 

Z  6i(r1  -  1) 

i-1 

s  ^  (l-«i)(yi.  -  yt)2 

S2(0)  -  ilUli _ _ _ L  .  (5) 

Z  (l-6.Hr  -  l) 

i-1 

It  can  be  seen  that  e(s2(1)  -  S2(0))  -  o2(l)  -  o2(0).  Thus  S2(l)  - 

2  2  2 

SjfO)  is  an  unbiased  estimator  of  a  (1)  -  a  (0). 

It  is  important  to  estimate  the  dispersion  effect  which  is  a  simple 

contrast  of  dispersions  at  levels  1  and  0,  or  in  other  words,  to  estimate 
2  2 

°  (tl»***»tu-l»  l»  CU+1 . O  "  0  (fcl . tu-l»  °*  tu+l»',,’tnJ  for  a 

(t i»*»*»tu_£*tu+jf ...,tu).  This  i».  possible  if  both  runs  (tj,...t  j, 

1»  tu+1 . tm)  and  [tl . cu-l»  °»  tu+l,***»tnJ  aPPear  as  rows  tn  T* 

Suppose  that  both  runs  are  present  in  T  and  they  appear  in  the  rows  1 . 
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2 

i»  of  T.  The  unbiased  estimator  of  the  parameter  of  interest  a  (_t_|  )  ~ 


o2(t!  )  is 
12 


i  i 

^(yi  i  -  y*  )2  ^2(y±  j  -  yA 

j-1  V  X1  _  j-1  2J  / 


^  "  *■ 
2 


Suppose  that  the  runs  (tj  1,  t^+1 , • • •  .t^ )  and  (c i » • • • »Cu-l * 

0,  t^+j,«.«,t')  also  appear  in  the  rows  i^  and  i^  of  T.  We  then  un- 

biasedly  estimate  the  dispersion  effect  o  (_t)  )  -  o  (t?  )•  The  linear 

3  4 

functions  of  the  dispersion  effects  —  [°2(t,'  )  _  o2(t!  )] 

2 

+  ~  [o2(l[  )  -  o2[t±  )]  and  1  [a2(t;  )  -  o2(tj  )]  -  ±  [o2(t|  )  - 
2  ^3  4  2  1  2  2  3 

a  (t!  )]  are  very  important  in  terms  of  the  "main  effects"  and  the 
4 

"interactions"  considered  in  Box  and  Meyer  (1986),  Srivastava  (1987), 
Note  that  the  first  linear  function  is  in  the  form  (3)  and  is  thus  again 
a  dispersion  effect  of  the  factor  chosen.  The  second  linear  function  is 
not  of  the  form  (3)  but  a  contrast  of  two  dispersion  effects  of  the  form 
(3),  There  is  a  possibility  for  an  array  T  that  the  unbiased  estimation 

of  o2(t’ ,  )  -  o2(t'f  )  is  not  at  all  possible  for  any  (t.,...,t  ,, 

12  2  2 
t  .,,«..,t  )  but  the  unbiased  estimation  of  o  (1)  -  a  (0)  is  still  pos- 

sible.  This  however  is  not  a  desirable  situation  and  the  corresponding 

ing  array  T  is  not  well  chosen.  In  section  3  we  discuss  this  issue  in 

detail. 


We  denote  ^  as  the  vector  of  observations  corresponding  to  the  runs 

with  the  chosen  factor  at  level  u,  u  -  1,0.  Notice  that  consists  of 

and  2s)'  Let  \  be  the  submatrix  of  X  corresponding  to  ^ be  the 

submatrix  of  £  corresponding  to  ^  be  the  ordinary  least  squares 

*  -1 

fitted  values  for  We  have  Xu(XrX)  X'^.,  u  •  1,0.  We  denote 


y  -  y  m  r  y  Rank  r  ■  V  and  V  S^(u)  -  v'r’r  y,  u  *  1,0.  Thus  V  S  (u) 
Ai^-utr-  uu  u  2  ^  u  ir  u2 

Is  the  sum  of  squares  of  the  ordinary  least  squares  residuals  at  level  u 
of  the  chosen  factor. 

Theorem  1.  We  have  for  u  *  1,0, 


a.  E(vuS2(u))  »  Trace  (^r^r^  +  V(l-u)Z(l-u)r(l-u)ri)’ 

b.  e(vuS2(u))  *  Trace  ruruEururu»  if  and  only  if  ^r^^j  *  0,  i.e., 

rlr*  -  0. 

Proof .  Since  we  have  assumed  no  significant  lack  of  fit  for  the  fitted 

model,  we  get  for  u  -  1,0,  Eli'r^r^)  »  Trace  r^ruE  «  Trace  ryE  r^.  It 

can  be  be  seen  that  the  submatrix  of  r^  corresponding  to  ^  and 

are  r  r'  and  r  rl,  . .  Thus  we  have  r  E  r'  =>  r  r'  Err' 
u  u  u  (1-u)  u  u  u  u  u  u  u 

+  r  r ,  VE,,  .r,,  »r'.  Hence  the  part  a  of  Theorem  1  is  true.  Since 

u  (1-u)  (1-u)  (1-u)  u 

Trace  r  r'  .E.,  .r,,  xr'  is  the  sum  of  squares  of  all  elements  in 

u  (1-u)  (1-u)  (1-u)  u 

1/2  1/2 

r  sE,,  s  and  furthermore  E.,  .  is  a  diagonal  matrix  with  all 

u  (1-u)  (1-u)  (1-u)  6 

diagonal  elements  positive,  the  part  b  of  Theorem  1  is  true.  This 
completes  the  proof  of  the  theorem. 

It  follows  from  Theorem  1  that  in  case  CjCq  »  0,  the  unbiased  estimator 


of  the  dispersion  effect 


Trace  r^'E^r*  Trace  V’E^r’ 


Trace  r^rj 


Trace  rQr' 


(7) 


is  ( Vj S2( 1)/Trace  r ^ r J )  -  (VQS2(0)/Trace  ^r^).  The  dispersion  effect 
(7)  is  of  the  form  (3). 

When  TjTq  *  0,  two  vectors  of  residuals  rg_  and  r^  at  levels  1  and 


0  of  the  factor  are  correlated  under  the  model  (1-2).  We  present  a 
vector  r(i-u)aZ.0^  "adjusted  residuals"  at  level  (1-u)  of  the  factor, 
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adjusted  w.r.C.  r^  (a  *  1,0)  so  that  rgay  and  r^y  are  uncorrelated 
under  (1-2).  Let  rui(vuxf0  be  a  submatrix  of  ru  so  that  Rank  rul  »  Vu, 


u  *  1,0.  We  write  for  u  *  1,0, 


r(l-u)a  =  r(l-u)^1  "  rul^rulruJ  rul^*  (8) 

It  can  be  seen  that  ruir(i_u)a  “  0  and  hence  rur(i_u)a  =  0*  can  be 

checked  that  Rank  r( i-u)a  a  ((N-p)  “  )  =  v(i_u)a  (say).  We  have  for 

u  =  1,0, 

r(l-u)a^  *  r(l-u)*-1  “  rul^rulrul^  rul}  r(l-u)-^(l-u)*  (9) 

Thus  for  u  =*  1,0,  T^i_u^aL  depends  on  only  through  and»  more¬ 

over,  Cov(rQa^,  r^^J  ”  0.  i*e*>  they  are  uncorrelated  under  (1-2).  Let 

r/t  ,  ,  be  a  submatrix  of  r,.  ,(l  -  r'.fr  ,r’  )  *r  ,lr!,  .  with  rank 

(l-u)al  ( 1 — u ) ^  ull  ul  ul'  ul'  (1-u) 

V(i_u)a>  u  “  1.0.  We  now  have  the  sum  of  squares  of  the  set  of  linear 
functions  r  [see  Scheffe  1959]  as 

SS(rual^  “  AlrualfrualruaJ'lrual^'  ( 10) 

with  d.f.  V  ,  u  =  1,0.  We  define  for  u  =*  1,0, 
ua 

It  can  be  seen  that  for  u  =  1,0, 

E!\jaS3a<“>]  ‘  ^(Wual^'ualWar  <12> 

Thus  the  unbiased  estimator  of  the  dispersion  effect 

Trace  (rlalrlal)~lrlalElrlal  Trace  ^r0alr0al r0al E0r0a 1 


is  (1)  -  (0).  Again,  Che  dispersion  effect  (13)  is  of  the  form 

3a  3a 


3.  Characterization  of  Arrays. 

In  this  section  we  observe  that  for  some  arrays  it  may  not  be 
possible  to  estimate  unbiasedly  both  o2( t^ p . . . , t^u  *  l,...,t^m)  and 

°  (til»*,,»tia  “  °*,***tiiJ  for  any  (til»***»ti(u-l)»ti(u+i) . tim^‘ 

We  present  the  characterizations  of  the  arrays  so  that  the  unbiased 

estimation  of  both  o2(t^p . . . ,t^u  =  l,...,t^  )  and 

o2( t ......  ,t ,  =  0,...,t,  )  are  possible  for  at  least  one 

v  ll  iu  im 

Ccil » • **ci(u— 1) » ci(u+l) ’ *  *  * ’ CinJ  and  for  every  u  ln  I1 . 


3.1.  Orthogonal  Arrays. 

We  first  consider  the  orthogonal  array  T  (nxm)  whose  rows 


_t ! ,  i  =  l,...,n,  satisfy 


At^  »  over  GF(2), 


where  A  (sxm)  is  a  known  matrix  of  rank  s  and  with  two  distinct  elements 


0  and  1,  c_  (sxl)  is  a  known  vector  with  at  most  two  distinct  elements  0 
and  1,  GF(2)  is  the  Galois  Field  with  two  elements  0  and  1.  We  have 
n  =  2ra  s.  Such  an  orthogonal  array  is  called  a  type-A  orthogonal  array 


(see  Srivastava  and  Chopra  1973).  We  write  the  matrix  A  in  (14)  as  A  = 

[a^l  ,a^2*  •  •  *  .a^] ,  where  a^sxl)  is  the  uth  (u  *  l,...,ra)  column  of  A. 

Theorem  2.  Consider  a  run  t'  «  (t.,...,t  ,  ,t  ,t  ,,,...,t  ]  of  the  type- 

-  —  v  1  u-lu’u+1  m 

A  orthogonal  array  T.  Then  o2(t, . t  . ,t  »  1,  t  t  1  and 

v  1  u-l  u  u+1  m 

a2(t. . t  ,,t  ”  0,  t  ,t  )  are  both  unbiasedly  estimable  if  and 

v  1  u-lu  u+1  m 

only  if  the  uth  column  of  A,  i.e.,  is  a  null  vector. 

Proof.  We  denote  t^  -  ( 1 1 , . . . ,  tu_L  ,tu  =  1,  tu+1 , . . .  ,tm)  and  t_'(0)  = 

^1 . Cu-rtu  *  °*  tu+l’*”*tm^  We  have  1(1)  +  1(0)  *  ^  over  GF(2)’ 

where  e^  is  the  vector  with  the  uth  element  unity  and  the  other  elements 

.  We  can  estimate  both  o2(t_^)  and  a2(t_^Q^)  if  and  only  if  both 


are  zero 
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and  appear  as  rows  in  T.  (We  know  that  one  of  them  is  already 

in  T.)  If  both  _t^  and  t_^  are  in  T  then  we  have  At^  =  At^Q^  -  c. 

We  get  At/1N  +  At/rtX  «  Ae  ,  i.e.,  c  +  c  =  a  .  But,  c  +  c  =  2c  -  0  over 

GF(2).  Thus  the  'only  if'  part  is  true.  The  'if'  part  is  trivial.  This 
completes  the  proof. 

We  want  to  estimate  a2(t_£^)  and  for  at  least  one  (t ^ ,.. . ,tu_j , 

t  +j,...,t  )  and  for  every  factor.  It  follows  from  Theorem  2  that  this 
is  possible  if  and  only  if  the  matrix  A  in  (14)  is  a  null  matrix.  Then 
the  type-A  orthogonal  array  T  does  not  exist.  We  therefore  conclude  that 

the  type-A  orthogonal  arrays  are  not  useful  in  the  unbiased  estimation  of 

the  dispersion  effects  of  the  type  o'  2(.L(i)]  °2U<o))  for  every  factor 

out  of  m  factors.  Although  the  type-A  orthogonal  arrays  are  widely  used 
in  quality  control  studies,  we  thus  observe  a  serious  drawback  of  such 
arrays. 

We  now  consider  the  orthogonal  array  T  (nxm)  whose  rows 
i  «  l,...,n,  are  solutions  of 

At_i  *  over  GF(2),  k  =  l,...,f  (15) 

where  A  (sxm)  and  c^  (sxl)  are  known  as  in  (14),  f  is  also  known.  We 
have  n  «  f  2m  s.  Such  an  orthogonal  array  is  called  a  type-B  orthogonal 
array  (see  Srivastava  and  Chopra  1973). 

Example  1.  Consider  a  25  factorial  experiment  with  8  runs  which  are  rows 


of  the  matrix  T  (8x5).  We  take 


X1  +  x2 


x3  +  x4 
x2  +  X4  +  x5 


Thus  s  -  3.  We  also  take 
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order  interactions  are  all  zero  (Resolution  III  or  Main  Effect  Plan). 

Theorem  3.  Consider  a  run  t'  *  (t,,...,t  ,  ,t  ,t  ,t  1  of  the  type 

-  —  v  1*  ’u-luu+l  m' 

B  orthogonal  array  T  satisfying  (15).  Then  a  (t,,...,t  »  l,...,t  )  * 

l  U  la 

o2(l(l))  and  02(t^f...,tu  *  0,...,tn)  *  °2(jL(o))  are  both  unbiasedly 

estimable  if  and  only  if  the  uth  (u  »  l,...,m)  column  of  A  in  (15)  is 

a  ■  c.  +  c.'  for  k,k'  in  ll,...,f}. 

— u  — k  — k  1  ’ 

Proof .  Me  have  t^  ^  +  1(0)  "  over  GF<2>*  111118  Ald)  +  Al(0)  ’  A^,  " 

a  .  Now  if  t...  and  t/ns  are  in  T,  we  take  At/IX  •  c,  and  At,„.  "  c. 

-u  -HD  — (0)  ^Kl)  -k  =-{l)  -k’ 

Hence  ^  ■  c^  +  c^*.  Suppose  now  that  a^  ■  c^  +  £k'"  have  to  s^ow 


that  both  t_£  and  t_^Qj  are  in  T. 
already  in  T  depending  whether  the 


We  know  that  one  of  and  t_(Q)  is 

uth  element  t  in  C_  is  equal  to  1  or 


4 


0.  Assuming  £^  ^  is  in  T,  we  show  that  ^.(q)  a^so  in  can  be 


seen  that 

A!(0)  “  A^u  "  1(1))  ”  iu  “  -^t  “  £k' * 

Thus  t^Q)  is  also  in  T.  The  rest  is  clear.  This  completes  that  proof. 

We  observe  that  if  a  type-B  orthogonal  array  satisfies  the  condition  of 

Theorem  3,  then  the  type-B  orthogonal  array  whose  runs  are  solutions  of 

At^  m  X.  ~  £k’  k  *  l,...,f,  i  =  l,...,n,  j_’  *  (1,1,.. .1),  also  satisfies 

the  condition  of  Theorem  3. 

Example  2.  We  again  consider  the  Example  1.  The  matrix  A  is  given 


below. 


1  1 


0  0  11 


1  1 


Note  that  none  of  the  columns  in  A  can  be  expressed  as  the  sum  of  £j  and 
£2«  Thus  for  the  type-B  orthogonal  array  T  in  (17),  the  condition  in 
Theorem  3  does  not  hold.  We  now  consider  a  2 5  factorial  experiment  with 
16  runs.  We  take  Ax  as  in  (16),  f  »  4,  £j  and  £2  as  i°  Example  1  and 


£3  -  0 


The  matrix  T  (16x5)  consists  of  all  runs  in  (17)  and  the  following  8 


additional  runs. 


/ 
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0  0  0  0 
0  0  0  0 
110  0 
110  0 
0  0  11 
0  0  11 
1111 
1111 


1 

0 

0 

1 

0 

1 

1 

0 


The  matrix  T  (16x5)  is  a  type-B  orthogonal  array  of  strength  1  but  not  of 
strength  2.  However,  it  is  nearly  an  orthogonal  array  of  strength  2. 

The  condition  of  Theorem  3  holds  for  T  (16x5).  Hence  °2(jL(i))  an<* 
a2(_t^o))  are  both  unbiasedly  estimable  for  a  (t ^  , .  • .  ,tu_^  »cu+p  •  •  •  ,tQ) 
and  for  every  chosen  factor.  We  thus  observe  that  we  need  more  runs  in 
the  array  to  satisfy  the  condition  of  Theorem  3  than  to  satisfy  the  con¬ 
dition  for  the  orthogonal  array  of  strength  2.  We  need  even  more  than  16 
runs  to  satisfy  both  the  condition  of  Theorem  3  and  the  condition  for  the 
orthogonal  array  of  strength  2. 

3.2  Balanced  Arrays  (B-arrays)  of  Full  Strength. 

We  now  consider  a  B-array  T  (nxm)  of  full  strength  with  n  distinct 
runs  (see  Srivastava  and  Chopra  1973).  We  first  describe  this  B-array. 
Let  be  the  set  of  all  (lxm)  vectors  having  i  elements  equal  to  1  and 

the  other  (m-i)  elements  equal  to  0,  i  -  0,1 . .  The  sets  , 

i  *  0,  l,...,m,  appear  as  rows  in  T  with  frequencies  X^ ,  where 

*  0  or  1,  i  ■  0,  l,...,m.  Clearly  the  number  of  vectors  in  is  [”) 


?  xt(p. 

i-0  1 


and  n 
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Theorem  4.  Consider  a  run  t^  »  ( c i » •  •  •  »tu-i  »cu  ,cu+l » *  *  *  ^  a  B-array 

of  full  strength.  Suppose  that  1  is  appearing  w  times  and  0  is  appearing 
(m-w)  times  in  t_' .  Then  a  (tj,...,tu  *  l*.*.»tffl)  *  o  f  1 )  ^  and 
o2(t1,...,tu  *  0,...,tffl)  «  °2(£(o)^  are  both  unbiasedly  estimable  if  and 
only  if  the  runs  in  S  and  S  .  .  where  x  *  1  if  t  *  0  and  x  »  -1  if 
t  =  1,  appear  as  rows  in  T. 

Proof.  Suppose  t  -  t/ix.  Then  the  runs  in  S  appear  as  rows  in  T.  It 
now  follows  from  the  B-array  of  full  strength  that  t_^g^  will  appear  as  a 
row  in  T  if  and  only  if  the  runs  in  appear  as  rows  in  T.  The  rest 

is  clear.  This  completes  the  proof. 

When  Sw  and  sw+x .  x  -  1  or  -1,  are  present  in  T,  we  can  estimate  un¬ 
biasedly  both  o2(t*,...,tJ_1,t*  -  ,t*)  and  o2(t^ . t*_p 

t*  =  0,t* ...... ,t*l  for  every  factor  u(u  *  l,...,m)  and  for  some 

u  u+1  ra 

(t|,...,tj_ptj+p...,tj).  If  the  runs  in  T  (nxm),  n  =  1  +  m,  are  the 

runs  in  Sq  and  S^  (or,  equivalently,  S^  and  Sm_j)  then  T  is  a  resolution 

III  plan  and  the  condition  of  Theorem  4  holds.  If  the  runs  in  T  (nxm), 

n  *  1  +  m  +  are  the  runs  in  Sq»  and  S2  (°r»  equivalently,  S^, 

S  ,  and  S  »)  then  T  is  a  resolution  V  plan  and  the  condition  of 
ra-i 

Theorem  4  holds  for  every  factor.  In  view  of  the  unbiased  estimation  of 

2  2 

the  dispersion  effect  o  -  o  B-arrays  of  full  strength  are 

very  good  having  smaller  number  of  runs.  The  resolution  IV  plan  where 
the  runs  in  T  (nxm),  n  -  2  +  2m,  are  the  runs  in  Sq,  S^,  Sj  and  Sm_j  is 
good  in  terms  of  estimating  both  location  and  dispersion  effects.  The 
resolution  V  plus  one  plan  where  the  runs  in  T  (nxm),  n  »  2  +  2m  +  (”), 


are  the  runs  in  Sq,  Sm,  S^,  Sm_j  and  S2  (see  Srivastava  and  Ghosh  (1976)) 
is  good  in  terms  of  estimating  both  location  and  dispersion  effects. 

4.  Final  Remarks 

In  the  model  (1-2)  and  throughout  the  paper,  we  present  the 

dispersions  in  the  general  form.  A  special  dispersion  model  assumes  the 
2  2.  2. 

additive  form  o  (tj,...,tm)  =  ( t j )+. . .+on (tQ ).  There  are  2m  dispersion 

2  2 

parameters  au(l)  and  ou(0),  u  =  l,...,m  under  the  additive  model.  The 
additive  form  is  implicitly  assumed  in  the  orthogonal  array  experimenta¬ 
tions  (see  Taguchi  and  Wu  1985).  The  additivity  in  the  orthogonal  array 
experimentations  simplifies  the  situations.  For  example,  we  have  for  the 

uth  (u  =  1 ,  • . .  ,m)  factor,  <j2(l)  -  o2(0)  -  o2(t_^jj)  -  a2(t^Qj)  -  cj2(1)  ~ 

2 

a  (0).  The  strengths  of  the  orthogonal  arrays  are  assumed  to  be  2  or 

more.  Under  the  additive  model,  the  unbiased  estimation  of  the  dis- 

2  2 

persion  effects  au(l)  -  au(0),  u  =  1 . .  is  always  possible  in  the 

orthogonal  array  experimentations. 


\ 


/ 
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